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1. Introduction 



For infrared (IR) divergent loop integrals, the integrand functions have non-integrable singularities 
through vanishing denominators. Based on dimensional regularization, the integral is expanded as 
a function of a parameter that approaches zero 

In this paper we report numerical results for the leading coefficients obtained by convergence 
acceleration (extrapolation) of a sequence of integral approximations. The methods are explained 
in detail in [Q], where also results are given using double precision arithmetic. In the present paper 
we give extensive results using quadruple precision and show that in many cases the accuracy can 
be improved to near the relative machine accuracy. 

According to the asymptotic behavior of the integral, we can explore linear or nonlinear ex- 
trapolation techniques. The asymptotic expansion gives rise to linear systems of the form 



n-l 

Sf;=Y,akq)k{£e), ^<£<n, (1.1) 
k=o 



where Se is generally a scaled version of the integral /(e), approximated numerically. A linear 
system solver or a linear extrapolation method can be used if the q)^ are known functions of eg. 
Otherwise, a nonlinear extrapolation may be suitable, depending on the nature of the (pt functions. 

For (small) £ > 0, the numerical integral approximation may be affected by singular integrand 
behavior which occurs at the boundaries and/or in the interior of the integration domain. The sam- 
ple integrals in this paper pertain to classes of one-loop vertex integrals which are two-dimensional, 
over the unit triangle {{x,y) \ <y < I —x < 1}. An iterated or repeated numerical integration 
can be performed efficiently with the general adaptive integration programs Dqags/Dqag from 



QUADPACK Q. 

The expansions derived symbolically in generally involve hypergeometric functions. We 
calculate the hypergeometric function numerically in Section jsj, using an extrapolation to handle 
the singularity in the integration interval. In Section ^ we present results for the case of one off- 
shell (p3 7^ 0) and two on-shell (p^ = p\ = Q) particles. The coefficients of the divergent terms 
in the integral expansion aie calculated with an extrapolation as the parameter e introduced by 
dimensional regularization goes to zero. The integrals in the sequence have integrand singularities 
on the boundaries of the integration region. 

Section ^ addresses IR divergent integrals with one on-shell and two off-shell particles, where 
integrand singularities may occur in the interior as well as on the boundaries of the integration 
domain. In this case, the integrals in the extrapolation sequence with respect to e involves an 
extrapolation to deal with the interior singularity. 



2. Extrapolation 

For a sequence {S{ei)}, which converges to the limit 5^ = Ymie^^QS{Ei), an extrapolation may be 
performed with the goal of creating sequences which convergence faster than the given sequence, 
based on its asymptotic expansion 

S(e) ~ ^ + Ai<pi (£) + A2<?)2(£) + • • . . (2.1) 
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Toi 
Til 



Tk--1,1 

TrtO Tk:_i^2 

ZkI 
Tic+1,0 



Tk-1 =0 
TjcO = Sk 

T|c,A+l = Tk-+1,A+1 + 



1 



Tk-+1,A — "^kX 



Table 1: e-algorithm table 

as £ ^ 0. In the context of series convergence we consider the limit of its partial sums. Some 
extrapolation methods allow summing divergent series to a value referred to as anti-limit. 

A linear extrapolation yields solutions to linear systems of the form 



S{ei) =ao + aiq)i{ef) + ...avq)v{ei), 



0,.. 



(2.2) 



of order (v + 1) x (v + 1) for increasing values of v g]. The sequence of ee may be geometric 
or another type of sequence that decreases to 0. As an example, Romberg integration relies on 
the Euler-Maclaurin expansion of the integral as a function of the step size e = h. Then is 
assumed to be an expansion in even powers of h, for the composite trapezoidal rule values S{h) with 
h = £>0. Values for oq^ y are obtained for successive v by solving the (v + 1) x (v + 1) 



systems of ( [2.2[ ) implicitly using the Neville algorithm. 

More general sequences of e include the sequence by Bulirsch, of the form l/bg with bg = 
2,3,4,6,8,12,... . (consisting of powers of 2, alternating with 1.5x the preceding power of 2). 
The type of sequence selected influences the stability of the process, which was found more stable 
with the geometric sequence than with the harmonic sequence (with the Bulirsch seuence in be- 
tween) [|^]. On the other hand there is a trade-off with the computational expense of S{e), which 
may become prohibitive for fast decreasing e . For the computations in subsequent sections we use 
scaled versions of be, e.g., bi/ 16. 

If the functions of e in the asymptotic expansion ( p| ) are not known, a nonlinear extrapola- 
tion or convergence acceleration may be suitable f jlTj [I^ , plj , ||]. As an example of a nonlinear 



extrapolation method, the £-algorithm [ ]17[ ] implements the sequence-to-sequence transformation 
by [ [l^ ] recursively; and can be applied when the q) functions are of the form <Pjfe(e) = £'^*log^*(£), 
under some conditions on Vt and jik and if a geometric sequence is used for £. The actual form of 
the underlying £-dependency does not need to be specified. 

Table [T] gives the recurrence of the £ algorithm for a sequence Sk, K" = 0, 1, . . . and depicts the 
layout of the computations in a triangular table. 
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Parameters 


Real/ 


Extrapolated 


ESTIM. REL. 




Total Int. 


1 


m 


n 


Imag . 


T? C CT IT T <^/^'{^) 

Kh^ U LI -.X ' ' 


IJ.KK. rLf 


V 


1 iME [^S ) 


1 


1 


1 


Real 


-1.4533220287861469626056457257e-02 


7.09e-29 


22 


7.58e-02 








Imag. 


-1.5079644737231007544620688243e-01 


2.34e-28 


22 


9.70e-02 


1 


2 


3 


Real 


8.41776716878103930()5426755543e-02 


2.85e-28 


22 


1.26e-01 








Imag. 


-2.29022 104446695927()8392670265e-0 1 


1.61e-29 


22 


1.20e-01 


2 


1 


1 


Real 


1 .08766468837090905 178932142e-02 


7.26e-26 


21 


1.08e-01 








Imag. 


2.63893782901542632030862046e-02 


1.66e-27 


22 


1.30e-01 


2 


3 


4 


Real 


-2.890568082315061 17361219566e-02 


2.38e-26 


22 


1.47e-01 








Imag. 


5.57897846432151278376445447e-02 


1.44e-27 


22 


1.42e-01 


3 


1 


2 


Real 


- 1 .02672 1 79830 1 7027333 1 993276e-02 


2.94e-27 


21 


2.42e-01 








Imag. 


-5 .65486677646 1 62782923279e-02 


4.29e-24 


21 


2.70e-01 


3 


4 


5 


Real 


8.12135882438810010895282e-03 


2.58e-24 


21 


1.79e-01 








Imag. 


- 1 .27463626362435298394 1 22e-02 


3.92e-25 


23 


2.13e-01 



Table 2: Integration and extrapolation results for 2F1 {l + l,l + mJ + m + n + l,z + iO) for relative integration 
error tolerance of 10^^^ and the Bulirsch sequence for extrapolation 



3. Hypergeometric function 

A representation of the hypergeometric function is given by the Gauss series 

°° r{a+k)r{b+k) 



F{a,b,c;z) = ^ 



k=Q 



r{c + k) k\ 



which has [ z [ = 1 as its circle of convergence and has an analytic continuation defined by the Euler 
integral 

Tic) t''-Ui-tY^''^^ 
= rwf5b)/„ (1-,,).' 

Rec > Reb > 0, which denotes a one- valued analytic function in the complex plane cut along the 
real axis from 1 to 00. 

For a numerical computation where z € R, we replace z by z + i5 and evaluate the limit of 



2F1 {a,b,c;z + iS) as 5 —> by solving linear systems of the form (2^) with (ptiSt;) = d^, 
1, . . . , V and S{de) = 2F1 {a,b,c;z + i5f), i = 0,...,v. 



Table 2 lists results for a problem set from [ ]10| ] where a = I + I, b = l + m, c = l + m + n+\ 
at z = 10 (real). The relative error tolerances for the outer and inner integration are set at 10^^^ and 
10^^^, respectively. Since b and c are positive integer, the numerator in the integrand is polynomial, 
so there are no end-point singularities. We use the general adaptive integrator Dqag of QuADPACK 
for the iterated integrations, with the 7-point Gauss and 15-point Kronrod rules applied on each 
subinterval in its adaptive partitioning strategy. 

Note that the weights and abscissae of the integration rules in DQAG are given to 33 digits 
and the relative machine accuracy in quadruple precision is of about the same order. Table 2 gives 
quadruple precision results, extending the (10-digit) accuracy of the double precision calculations 
reported in [^. As compiler we use the intel Fortran Composer XE with the -rl6 flag. The cal- 
culations are run on a Macbook-Pro laptop with 3.06 GHz Intel Core 2 Duo processor and 8 GB 
memory (Mac OS X Version 10.6.4). 



According to (2.2) we obtain =5^^^' « c?o for the systems of order (v+l)x(v + l), v = l,2, 



The difference of successive results in £"1^^ = |^(^) — gives a measure of convergence. 
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Table 2 lists l,m,n followed by the result the relative measure of convergence Er = 

Eci^^ /\^^^^\, the value of V and the time taken (in seconds) for the integrals S{5e), i = 0,. . . ,v, 
needed to obtain the result. The other times in the program are not taken into account, including the 
system solving time which was measured and found negligible compared to the integration times. 

While we could have listed the result obtained after the error falls below 10^^^, we instead 
listed the result with the smallest estimated relative error. As expected, better accuracy is obtained 
for the smaller values of l,m,n. 



4. Asymptotics for one off-shell, two on-shell particles 

We first address the IR divergent integral ^3 {p\,P2,pl',njf,ny) from with one off-shell (p^ / 
0) and two on-shell particles {pj = p^ = 0), which we denote here by jT,{p\,p\,pl;n^,ny;e) = 



eT{-e) 



Anjxl J -pi nx + ny + 2e 



where pa^ = P3 + iO and is a renormalization constant (which we will replace by /x| ^ / {4n) , 
Ye is Euler's constant). The introduction of the parameter £ pertains to the dimension regularization 



technique [13] 



When = ^7 7^ or «y / 0, we have 



I^^\e) - \ +Co + ^{e) ] with (4.2) 



P3 



e 



1 

C I — — , Co 




A linear extrapolation can be formulated using systems of the form ( |2.2| ) with Se{ee) = £el{se) 
where I{e() « /3'^'*'(e) and (pk{ee) = e^. 

Table ^ displays results for the real parts of the coefficients C i and Co in the expansion, for 
nx = Tj =2 and riy = 0. In view of the integrand singularity along y = through the factor y^^^ 
where the exponent approaches - 1 , and a second derivative singularity along x = through the fac- 
tor in (p~Tl), we perform the calculation of I{ei) with the program Dqags of QuADPACK [|l^], 
which is equipped to deal with these end-point singularities. DQAGS uses the 7-point Gauss and 
15-point Kronrod rule on each subinterval created in its adaptive subdivision strategy. 

We set the requested relative accuracies to 5 x 10^^^ for the outer and 10^^^ for the inner 
integral. Table ^ lists the sequence of extrapolation results and shows how much accuracy can be 
reached. (After that, the accuracy no longer improves; i.e., stagnates for a few steps and decreases). 
The difference between successive results provides a good estimate of the convergence. For this 
problem it usually gives a somewhat conservative bound for the actual error, in the sense that the 
difference with the result of the previous step is in fact a measure of the error in the previous step. 
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1 

2 
3 
4 
5 
6 
7 
8 
9 
10 
11 
12 
13 
14 
15 



Integration 



# EVALS 



1125075 
1349535 
1348875 
134 9355 
1324125 
1274625 
127 4625 
1225125 
1200375 
1175625 
1101375 
1076625 
1037475 
1037475 
952425 



Total time : 



Time 



1.10 
1.34 
1.33 
1.33 
1.30 
1.26 
1.25 
1.21 
1.18 
1.16 
1.08 
1.06 
1.02 
1.00 
0.90 



20.31 



Extrapolation Results 



0. 8707294 1 267927795009280 1 1 356924 
0.3457684478555631741492957488633 
0.5262024637357823738306365020734 
0.4976817315764365660976698223580 
0.5001453140694736901417736134445 
0.500000 1 92 1 1 3 1 322347567 1 04032762 
0.49999999594673227 1 3652940448364 
0.4999999999993260163866288358759 
0.5000000000000053069685 193 1745 12 
0.4999999999999999722450728860463 
0.500000000000000000 1 08902333654 1 
0.4999999999999999999997 1 5 88863 1 9 
0.5000000000000000000000005564493 
0.4999999999999999999999999992988 
0.4999999999999999999999999999980 



Exact : 



0.5 



Co 



-0.5511174754172542062710316596 
3.1236092783487492253335060482 
0.77796707 1 9058996294760762565 
1.3769024472521615918683765304 
1.3030216329875385877494858112 
1.3025651116288535549087530061 
1 .302585079 1 32 1435704903401609 
1 .3025850932043285642085 160900 
1 .302585092994063225 1 1 667 1 6065 
1 .302585092994063225 1 1 667 1 6065 
1.3025850929940455873126312708 
1 .3025850929940456843793808956 
1.3025850929940456840169987640 
1.3025850929940456840179932776 
1.3025850929940456840179914927 



1 .3025850929940456840 1 799 1 4547 



Table 3: Integration performance (Dqags) , for rel. integration error tolerances 5x10 (outer), 10 
(inner); n^ — T] =2, ny = and = 100. Extrapolated real values. 



^-25 



The final accuracy reached for C i is about 2 x 10^ (absolute error), 4 x 10^ (relative error). 
The accuacy in Co lags behind and reaches about 4 x 10^^^ (absolute) or 3 x 10^^^ (relative error). 
When rijf = riy = we have 

^3°'°(^)-g(|? + %^+C0+^(£)) With 

C_2 = l, C_i=ln(-;;2)^ Co = -^ + Uvl\-pI). (4.3) 

With relative tolerances of 10^^^ for the outer and 5 x 10^^^ for the inner integrations, we obtain 
the best relative accuracies for C 2, C-\ and Co at v = 15, of 7.88e-25, 4.33e-22 and 4.0e-19, 
respectively. Note that the best accuracies for this problem obtained in double precision with 
integration tolerances of 10^'^ (outer), 5 x 10^^^ (inner) are reported in as 3.06e-12, 1.44e-10 
and 7.99e-09 at v = 11. The latter can also be compared with our 10^^^, 5 x 10^^^ quadruple 
precision result at v = 1 1, which yields 3.18e-14, 3.04e-12 and 5.24e-10, respectively. 

5. Asymptotics for one on-shell, two off-shell particles 

In the case of one on-shell {p\ = 0) and two off-shell {p\ / 0, 773 / 0) particles, the IR diver- 
gent integral is 

h{0,plpl,n,,n,;£) = -^h{0,0,p\,n,,ny-e) 2F1 (1, 1 - £,2 + n,; ^^^) (5.1) 
(471)-^ ^ p-i rix+l 

_ 1 £r(-e) fi-^ x"'/>' 

"(W(4^yo ""L \-{pl-pl)xy-ply{\-y)-iOy->^- 

(5.2) 
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For = Uy = we use the expansion 



73 {0,plpl;n^,ny;a) ~ +Co + ^(e) (5.3) 



with 



^ 1 In(l-z) , In^(l-z) 



2_ 2 

with z = '^^pj^ ^nd the C i coefficient of In the subsections below we will give results for 
two methods for the computation of ^3 {0,P2,p^',n_x,ny;e), the first based on (5J) and the second 



on 

5.1 Computation of /3 (0,pl,pl;nx,ny) with hypergeometric function 

We intend to perform an extrapolation with respect to the parameter e for dimensional regu- 



larization. According to (5. 1 ), each term in the extrapolation sequence consists of a double integral 
multiplied with a hypergeometric function, both dependent on e. We calculate the hypergeometric 
function numerically using the 5 extrapolation of Section ||, to alleviate the characteristic non- 
integrable singularity on the real axis when it occurs inside the integration interval. 



^_ 2 

Here z = f ^ in the fourth argument of 2F1 in ( 5. 1 ) is replaced by z + id . Thus in this process 



we need to evaluate sequences of hypergeometric functions of the form 

2 2 

2^1 (l,l-£^,2 + n,; ^^~^^ +i5^), jc=l,2,... 

for e = 1,2, The number of K"- values needed depends on the convergence for each (fixed) e. 

As the exponents b —\ = —Ee^ and c — b—\=nx + £i in the integrand numerator f'^^' (1 — tY^'^^^ 



of the Euler integral (3.1) are non-integer, we use the QUADPACK program Dqags to treat the 



integrand behavior at the end-points. 

Extrapolated (real part) results for C_i and Co in (5.3) and actual absolute errors are listed in 
Table 0. By setting p\ = 40, p\ = — 100, for the purpose of a numerical example, we have that 
Re{z) = 1.4, so that the hypergeometric function has an integrand singularity at ? = 1/1.4 in the 
interior of the integration interval (0, 1). 

The relative error tolerances for DQAGS were 10^^^ for the outer integration and 5 x 10^^^ for 
the inner integrals. The maximum number of subdivisions was set to 150 for 73 in both directions 
and to 300 for 27^1. For the dimensional regularization extrapolation, the Bulirsch sequence was 
used starting at 3, i.e., 3,4,6, 



For p| = —100, the integral in (11) is real and is multiplied with the (complex) value of 
the hypergeometric function. The real and imaginary parts of the latter are currently computed 
separately which, for the real part, took between about 0.44 and 0.55 seconds, and for the imaginary 
part between 0.33 and 0.41 seconds. Thus the time for the hypergeometric function computation 
for each equation in the linear system is under a second and is dominated by the times for the 
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Time 




Extrapolation Results 




V 


INT.(S) 






IErrorI 


Co 


IErrorI 


4 


3.47 


-4.2192812666950192419664334e-05 


7.47e-07 


-3.7121279333017303070237e-04 


2.39e-05 


5 


3.47 


-4 


1 395 8894045400495 11801 909e-05 


5.04e-08 


-3 .975 1 1 26098970774666 1 84e-04 


2.42e-06 


6 


3.41 


-4 


1448485840565542759037853e-05 


2.21e-09 


-3.9493403562445857754728e-04 


1.59e-07 


7 


3.39 


-4 


1446205836959577424975383e-05 


7.16e-ll 


-3.9510047588769404693383e-04 


7.45e-09 


8 


3.33 


-4 


1446278988 1 82334267440660e-05 


1.53e-12 


-3.9509279500930457847498e-04 


2.32e-10 


9 


3.30 


-4 


1 446277437326639980970560e-05 


2.43e-14 


-3 .95093032290225 8043049 1 e-04 


5.25e-12 


10 


3.53 


-4 


1446277461859729880198253e-05 


2.56e-16 


-3.9509302696654529617250e-04 


7.97e-14 


11 


3.86 


-4 


1 44627746 1 602 1 577 1 7453434e-05 


2.01e-18 


-3.95093027047165383 11163e-04 


8.86e-16 


12 


3.75 


-4 


1 44627746 1 604 1 8243208 1 1 46e-05 


1.05e-20 


-3.9509302704627248396080e-04 


6.66e-18 


13 


3.81 


-4 


1 44627746 1 604 1 7 1 84995795 2e-05 


4.14e-23 


-3.9509302704627918244479e-04 


3.67e-20 


14 


3.77 


-4 


1 44627746 1 604 1 7 1 89 1 430927e-05 


1.08e-25 


-3.9509302704627914557531e-04 


1.37e-22 


15 


3.81 


-4 


1 44627746 1 604 1 7 1 89 1 3225 14e-05 


3.19e-28 


-3 .95093027046279 1 457 1 332e-04 


5.44e-25 


16 


3.76 


-4 


1 44627746 1 604 1 7 1 89 1 32288 1 e-05 


4.87e-29 


-3 .95093027046279 1 457 1 267e-04 


1.12e-25 


17 


3.82 


-4 


1 44627746 1 604 1 7 1 89 1 322692e-05 


1.39e-28 


-3 .95093027046279 1 457 1 3 1 5e-04 


3.69e-25 


18 


3.78 


-4 


144627746 1 604 1 7 1 89 1 323096e-05 


2.63e-28 


-3.9509302704627914571 171e-04 


1.07e-24 


19 


3.82 


-4 


144627746 1 604 1 7 1 89 1 322859e-05 


2.63e-29 


-3 .95093027046279 1 457 1 292e-04 


1.40e-25 


fie ; 


-4. 1 44627746 1 604 1 7 1 89 1 322832e-05 




-3 .95093027046279 1 457 1 278e-04 





Table 4: Integration performance (Dqags)^, for rel. integration error tolerances 10^^^ (outer), 5 x 10^^^ 
(inner); vertex with one on-shell, two off-shell particles, nx — ny — and P2 — 40, — —100. Extrapolated 
real values and corresponding actual absolute errors. 




pmditced by GRACEFIG 



Figure \: qq ^ qqg diagram. The vertical gluon propagator carries p\ and the u-quark emitting a gluon 
carries p\. 

integrations of Jt, in ( pTT| ) (given in Table The total time spent in the double integration is 68.3 
seconds (including that of the first three steps, which are not shown). 

A sample diagram for the vertex correction with one on-shell and two off-shell particles is 
depicted in Figure |l| for the in qq qqg interaction. Here the vertical gluon propagator carries p\ 
and the u-quark emitting a gluon carries p\. Note that the vertical gluon propagator can be virtual 
so that the square of its momentum becomes negative. 



5.2 Computation of {0,pl,pl;nx,ny) using direct integration 



In this section we calculate the integral in (5.2) directly, which becomes highly singular at 



J = as £ — > and non-integrable for £ = 0, in view of the factor in the integrand 



M5F 



It 
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can be seen that {pj — p\)x + p2{'\-—y)=0 along a line that goes through {P2/ {p\ — pi) ,Q) and 
(0, 1). Thus for pI = 40, p\ = -100, this line runs through the integration domain, where it causes 
a singularity that becomes non-integrable as £ ^ 0. 

According with previous work (e.g., ||^ 0, ^, 18, 3|]), we will replace /O by i5 in the integrand 



denominator of (5.2). This leads to integrals of the form 



/3(£,5) 



1 £r(-£) 



(47r)2 [An^l 



2\£ 



dx 



dy 



i-ipi 



-pl)xy-ply{l-y)-i5)^ 



+ Co(5) + ^(£) 



(5.4) 



to be computed in an extrapolation as 5 ^ 0, for each (fixed) value of the dimensional regular- 



ization parameter £. The expansion ( p. 3D or ( |5.4| ) with respect to £ justifies a linear extrapolation. 
However, a linear extrapolation as 5 —> is not suited, as it cannot be assumed that the underlying 
expansion is in integer powers of 5. 

Preliminary results are shown in Tables for the real and the imaginary parts, respectively. 
We employ a sequence of 5^ = 2^^^'^, K" = 0, 1, . . . in an extrapolation using the £-algorithm of 
Wynn [15, 17]. The current implementation incorporates a simple strategy with an extrapolation 
sequence of fixed length (= 18), which should be improved with a termination criterion based on 
the estimated errors. Geometric sequences of 5 with a smaller ratio should also be tested. For the 
extrapolation with respect to £ the Bulirsch sequence 3,4, 6, . . . is used. 



V 


C-i 


IErrorI 


Co 


IErrorI 


4 


-4.2192812666986e-05 


7.47e-07 


-3.712127933292e-04 


2.39e-05 


5 


-4.1395889404493e-05 


5.04e-08 


-3.9751 12609915e-04 


2.42e-06 


6 


-4.144848584()595e-05 


2.21e-09 


-3.949340356225e-04 


1.59e-07 


7 


-4.1446205836943e-05 


7.16e-ll 


-3.951004758891e-()4 


7.45e-09 


8 


-4.1446278988175e-05 


1.53e-12 


-3.950927950097e-04 


2.23e-10 


9 


-4.1446277437321e-05 


2.43e-14 


-3.950930322905e-()4 


5.24e-12 


10 


-4.1446277461827e-05 


2.23e-16 


-3.950930269725e-()4 


7.38e-14 


11 


-4.1 44627746 1600e-05 


3.96e-18 


-3.950930270436e-04 


2.69e-16 


Ex: 


-4.1446277461604e-05 




-3.95093027()463e-()4 





Table 5: Integration performance (Dqags)^, for rel. integration error tolerances 10^'^ (outer), 10^'^ 
(inner); Vertex with one on-shell, two off-shell particles, n.^ — ny — and P2 — 40, = — 100. Extrapolated 
real parts and actual absolute errors. 



V 


C-i 


IErrorI 


Co 


IErrorI 


1 


0.758010428109e-04 


6.63e-05 


9.62356140022e-04 


4.38e-04 


2 


1.464854199104e-04 


4.38e-06 


4.67565500325e-04 


5.66e-05 


3 


1.420636477897e-()4 


3.90e-08 


5.25048537895e-04 


4.91e-07 


4 


1.421 00046009 le-04 


2.58e-09 


5.24284175287e-04 


8.47e-08 


5 


1.421026578481e-04 


3.01e-ll 


5.24197984601e-04 


1.48e-09 


6 


1.421026279735e-04 


2.13e-13 


5.24199448458e-04 


1.55e-ll 


7 


1.421026277585e-04 


2.16e-15 


5.24199464154e-04 


2.32e-13 


8 


1.421 0262776 12e-04 


5.84e-16 


5.24199463865e-04 


5.64e-14 


Ex: 


1.421026277606e-05 




5.24199463922e-04 





Table 6: Integration performance (Dqags)^, for rel. integration error tolerances 10^'^ (outer), 10^'^ 
(inner); Vertex with one on-shell, two off-shell particles, n^^ = ny = and p\ = 40, p\ = — 100. Extrapolated 
imaginary parts and actual absolute errors. 
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The requested accuracies for Dqags are 10~ for the outer and 10~ for the inner integra- 
tions. The real and imaginary parts are calculated separately as the QuADPACK programs currently 
do not handle complex functions. The imaginary parts converge more quickly than the real parts 
for this problem. We list the steps 1 to 8; the last for a 9 x 9 system. The results do not improve 
after that. Note that it is not necessary to solve the previous systems, but the integrals need to be 
calculated to set up the system. Successive systems can be solved to provide an estimate of the 
error. 

Conclusions 

We use ID integration programs from QUADPACK for the iterated integrations underlying the ex- 
trapolation processes for IR divergent vertex integrals. We give preliminary results in quadruple 
precision which demonstrate that the calculations result in high accuracy. For the vertex computa- 
tion with one on-shell and two off-shell particles we devise a double extrapolation, since problems 
are introduced with integrand singularities in the interior of the domain as well as on the bound- 
aries. Many integrals result, which makes the computation a good candidate for parallelization, 
especially when applied to more complicated diagrams and possibly in higher precision. 

More work is needed for improvements and extensions of these strategies and analysis of 
the numerical methods and results. This computation is a step toward a more automatic numerical 
handling of various types of loop integrals, thereby circumventing the need for a precise knowledge 
of the location and structure of the singularities. 
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